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Mixing 
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An overview of selected experimental results in the field of D^-D^ oscillations is presented. The average results 
for the mixine parameters, x = (0.89 ± 027)% and y = (0.75 ± o{g)%, exclude the no-mixing hypothesis at 
the level of 6.7 standard deviations. No sign of CP violation in the -D" system is observed. The measurements 
impose constraints on the parameter space of many New Physics models. 
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1. Introduction 

Last year, 31 years after the discovery of D° mesons, 
the first evidence of a mixing phenomena in the sys- 
tem of neutral charm mesons has been obtained [l|, Q ■ 
Following this breakthrough were several additional 
measurements enabling - through an averaging proce- 
dure - a quite precise determination of the parame- 
ters governing the mixing. The results presented in 
the paper follow from the data collected by the two 
B-factories experiments, Belle and BaBar, from the 
charm-factory experiment Cleo-c, as well as from the 
proton collider experiment CDF. At the B-factories 
the cross-section for the continuum production of cc 
pairs is around 1.3 nb, which with the integrated lu- 
minosity of KEKB amounts to 10^ produced charmed 
hadron pairs. At CESR, the --800 pb^^ data sample 
corresponds to 2.8 x 10® D^D^ pairs produced in a co- 
herent C = — 1 state. And while at the Tevatron the 
experimental environment for the presented measure- 
ments is more difficult, the cross-section for neutral 
charm meson production with a transverse momentum 
larger than 5.5 GeV/c yields a starting data sample 
of 50 X 10^ D^'s. The experiments thus provide a re- 
ally diverse experimental environment for successful 
studies in charmed hadron physics. 

The mixing, that is the transition of a neutral 
meson into its antiparticle and vice-versa, appears as 
a consequence of states of definite flavour (Z?°, 
being a linear superposition of the mass eigenstates 
(states of simple exponential time evolution, D12)' 



\Di,: 



(1) 



It is governed by the lifetime of D mesons, r = l/F, 
and by the mixing parameters x = (mi — 77x2) /F and 
y = (Fi — F2)/2F. mi^2 and Fi_2 denote the masses 
and widths of the mass eigenstates Di and D2, re- 
spectively. F — (Fi -I- F2)/2 is the average decay 
width. The mixing rate is small. The contribution of 
loop diagrams, successfully describing the oscillations 
in other neutral meson systems, is suppressed since 
the D'^ system is the only neutral meson system with 
down-like quarks exchanged in the loop. In this short 
distance calculation x is negligibly small due to the 
small SU(3) flavor symmetry breaking (rnj. ~ rn^ ^) 



matrix {\Vub\ ~ 4 x 10"'^) [3|. Long distance contribu- 
tions to the D^-D^ transition are difficult to calculate. 
Current theoretical estimates predict the mixing pa- 
rameters \y\ < 10~^ [lilll- 

Violation of the CP asymmetry (CPV) in the 
charm sector is expected to be small. Since the pro- 
cesses with mesons involve mainly the first two 
generations of quarks, for which the CKM elements 
are almost real, the expected level of CPV is 0(10"^) 
which is below the current experimental sensitivity. 



2. Measurements 

Several methods and selection criteria are common 

to the presented measurements. Tagging of the flavour 

of an initially produced meson is achieved by 

reconstruction of decays D^^^t O'" ~^ 

— _ 

D TTg . The charge of the characteristic low momen- 
tum pion TTs determines the tag. The energy released 
in the D* decay, 



q = M(D*)- M{D^)- 



(2) 



and elements of Cabibbo-Kobayashi-Maskawa (CKM) the other hand, flavour speciflc final states like K 



has a narrow peak for the signal events and thus helps 
in rejecting the combinatorial background. Here, 
M{X) is used to denote the invariant mass of the X 
decay products, and mx stands for the nominal mass 
of X . mesons produced in B decays have different 
decay time distribution and kinematic properties than 
the mesons produced in fragmentation. In order to 
obtain a sample of neutral mesons with uniform prop- 
erties one selects D* mesons with momentum above 
the kinematic limit for B meson decays (B-factories) 
or uses the impact parameter distribution to isolate 
primary charm mesons (CDF). 

2.1. Decays to CP eigenstates 

In the limit of negligible CPV the mass eigenstates 
_Di.2 are also CP eigenstates. In decays fcp 
only the mass eigenstate component of 13° with the 
CP eigenvalue equal to the one of fcp contributes. 
By measuring the lifetime of D'^ in decays to fc p one 
thus determines the corresponding 1/Fi or I/F2. On 
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have a mixed CP symmetry. The measured value of 
the effective hfetime in these decays corresponds to a 
mixture of l/Fi and \/T2- The relation between the 
two lifetimes can be written as |g| 



rifcp) 



1 



VfVCP 



(3) 



where T{fcp) and t{D'^) arc the lifetimes measured in 
fcp and D'^ K^tt^, respectively, r/f — ±1 
denotes the CP eigenvalue of fcp- The relative dif- 
ference of the lifetimes is described by the parameter 
ycp- 

CP-even final states fcp = ,tt^tt^ were 

used to measure ycp [H- Expressed in terms of the 
mixing parameters, ycp reads @ 



1 



ycp = y cos ( 



-Am sin ( 



(4) 



with Am and (j) describing the CPV in mixing and 
interference between mixing and decays, respectively. 
In case of no CPV {Am,4> = 0) and ycp = U- 

Simultaneous fits to decay time distributions of se- 
lected 1?° K~ , K^TT'^ and tt+tt" candidates 
were performed with ycp as a common free param- 
eter. The fit is presented in Fig. [Tlla)-(c). The agree- 
ment of the fit function with the data is excellent, 
X^/n.d.f = 312/289. The final value obtained is 



ycp (1.31 ±0.32 ±0.25)% 



(5) 




Figure 1: (a)-(c): Result of the simultaneous fit to decay 
time distributions in D^' decays to KK, K-n and tttt final 
states. The hatched areas represent the contribution of 
backgrounds, (d): Ratio of D° ^ fcp and D" K~tv+ 
decay time distributions. The slope visualizes the difi'er- 
ence of effective lifetimes. 



The largest systematic uncertainties arise from the 
assumed resolution function (common offset in indi- 
vidual decay modes), possible deviations of accep- 
tance dependence on decay time from a constant (esti- 
mated by a fit to the generated t distribution of recon- 
structed MC events) and variation of selection criteria 
(effect estimated using high statistics MC samples). 

The resulting ycp is more than 3 standard devia- 
tions above zero and hence represents clear evidence 
of — D mixing, regardless of possible CPV. The 
difference of lifetimes is made visually observable by 
ploting the ratio of decay time distributions for decays 
to fcp and K~7r+ in Fig. md). 

Recently the BaBar collaboration performed a sim- 
ilar measurement 0, with results for individual life- 
times shown in Fig. [21 The obtained value of the 
mixing parameter is 



ycp ^ (1.24 ±0.39 ±0.13)% 



(6) 
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Figure 2: Lifetimes of D'^ mesons measured in decays to 
Ktv, KK and tttt (separately for D^' and D'^ for the latter 
two modes) 0]. The average value of the lifetime as mea- 
sured in decays to CP eigenstates is lower than the one 
measured in flavour specific final state. 

The two measurements make ycp the most pre- 
cisely measured individual mixing parameter in the 
D'^ system. 

Final states of definite CP allow also a search for 
possible CPV. Two methods, decay time dependent 
and time integrated, have been exploited. For the for- 
mer, the lifetime in D'^ — *■ /cf is measured separately 
for D'^ and tagged events. The asymmetry is @ 



Ar - 



t{D fcp)-T{D° ^ fcp) 



t{D 



-AmJ/COS 



fcp) + TiDO 
X sin <h 



CP) 



(7) 



The values of Ar measur ed by Belle ^ and BaBar 
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are 

At = (0.01 ±0.30 ±0.15)% 

At = (0.26 ±0.36 ±0.08)% , (8) 

respectively, and show no sign of CPV at the level of 
around 0.3%. 

With the time integrated method one measures the 
asymmetry 

r(i?" ^ fcp) - r(^o fcp) _ 



r(i?o ^ fcp) + r(^o ^ /cp) 



(9) 

ylcp receives contribution from all three types of 
CPV , direct, CPV in mixing and in the interference 
between decays with and without the mixing. The 
latter two are independent of the final state. Experi- 
mentaly the measured asymmetry must be corrected 
for possible charge asymmetries in the detection of 
the slow pion as well as the forward-backward asym- 
metry {Afb) in the production of fermion pairs in 
e'^e~ collisions. The method of determination of tt^ 
correction factors was developed in ^ using the un- 
tagged K~Tr~^ decays. The forward-backward 
asymmetry is separated on the basis of its symmetry 
properties as a function of the D meson polar angle in 
the center-of-mass system. Figure [3] shows the mea- 
sured Acp ((a),(b)) and Afb ((c), (d)) as a function 
of the D meson polar angle Averaging over the 
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Figure 3: Time integrated CP asymmetry ((a),(b)) and 
forward-backward asymmetry ((c),(d)) as a function of the 
D meson polar angle @]. 



polar angle yields the value 

= (0.00 ±0.34 ±0.13)% . (10) 
Measurement of the Belle collaboration j9|] 

= (-0.43 ±0.30 ±0.11)% . (11) 
is also consistent with no CPV. 



2.2. Wrong-sign decays to hadronic final 
states 

Decays of D mesons to two-body hadronic final 
states accessible to both, D° and Z)°, have tradi- 
tionally been used to search for the mixing. Final 
state K^TT^ can be reached through a doubly Cabibbo 
supressed (DCS) D'^ decay as well as through the 
D'^ — > Z)° mixing followed by a Cabibbo favoured 
(CF) decay. The time evolution for these decays has 
three terms: 

(12) 

The first one is due to DCS decays, the third one 
due to the mixing, and the middle term represents 
the interference of the two contributions. Rd is the 
Cabibbo suppression factor relative to CF decays. 
x' and y' are the mixing parameters, rotated by a 
strong phase difference between DCS and CF decays, 
x' = xcosS + ysinS and y' = ycosS — xsinS. The 
dimensionless time t is measured in units of t{D'^). 

The BaBar collaboration obtained the first evidence 
for D'^ mixing by performing the decay time study of 
wrong-charge decays D*^ D'^tt'^, K^tt^ 
to separate the DCS and the mixing contribution. 
The parameters obtained from the fit are presented 
in terms of likelihood contours in Fig. [4] (top). The 
central value lies slightly in the non-physical region 
{x'^ < 0) and the no-mixing point {x'^ = 0,y' = 0) 
is excluded at the level corresponding to 3.9 standard 
deviations. Recently CDF collaboration obtained the 
result of similar significance [l3| shown in Fig. [J] (bot- 
tom left). Result form the Belle collaboration [ll| 
takes into account the presence of a physical bound- 
ary and is presented in Fig. |4] (bottom right) as a 95% 
C.L. contour calculated using the Feldman-Cousins 
method. 



2.3. Time dependent Dalitz analyses 

Several intermediate resonances can contribute to 
a hadronic multi-body final state. In a specific de- 
cay channel D'^ — s- Kstt'^tt^, recently analyzed by 
Belle fl3], contributions from CF decays (e.g. 
K*^n'^), DCS decays (e.g. K*^n~) and decays 

to CP eigenstates (e.g. D'^ p'^Ks) are present. In- 
dividual contributions can be identified by analyzing 
the Dalitz distribution of the decay. Moreover, for 
a self-conjugated final state these different types of 
decays interfere and it is possible to determine their 
relative phases (unlike in the case of Z?" if+7r~ 
decays). Since these types of intermediate states also 
exhibit a specific time evolution one can determine di- 
rectly the mixing parameters x and y by studying the 
time evolution of the Dalitz distribution. 
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Figure 4; Top: 1 - 5 cr likelihood contours for parameters 
a;'^ and y' obtained from a fit to the decay time distribu- 
tion of -D° — > K^-K~ decays @]- Bottom left: Bayesian 
probability contours in x'^ , y' plane corresponding to 1 - 4 
a as obtained in [lo|. Bottom right: 95% C.L. contours of 
x''^,y' using Feldman-Cousins approach from pH ]. 



The signal p.d.f. for a simultaneous fit to the Dalitz 
and decay-time distribution is 

M{m't,m\,t) = {KsTT+7r-\D"{t)) = 



+ ^A{m^,,ml)[e 



(13) 



The matrix element is composed of an instantaneous 
amplitude for D'^ decay, A{m^,ni^), and an am- 
plitude for the D° decay, A{m'^,m-^), arising due 
to a possibility of mixing. They both depend on 

M^iKsn^) and m\ = 



the Dalitz variables 



no CPV (Stat, only 
no CPV 

CPV (Stat, only) 
CPV 




x(%) 



Figure 5: 95% C.L. region for parameters x and y as ob- 
tained in D'^ — > KsTr'^'K~ decays p^ . 



28 

6 = (-0.24 ± ' ±0.09)rad . (14) 
^ ^ 0.30 ' ^ ' 

The measurement represents the most accurate deter- 
mination of a;. The CP violating parameters \q/p\ and 
(j) = arg{q/p) are consistent with no CP violation. 

BaBar has performed the time dependent Dalitz 
analysis of D° K'^tt~tt° decays [l3|. The final 
state is flavour specific and hence the wrong-sign de- 
cays again receive contribution from mixing and DCS 
process. The combined Dalitz-decay-time signal dis- 
tribution has the form 

\{K+7T+n-\D"m'^[\Af\' + 
+ \Af\\Af\{y" cosSf - x" smSf)t + 



^t^\e-' 



(15) 



A J (depending on the Dalitz variables AP{K^tt ) 
and M^{K+7r°)) is the amplitude for K+tt~tt° 
decays determined from the fit to the Dalitz distri- 
bution of wrong-sign decays. The amplitude for Z)*^ 
decays, Ajr, is fixed to the values obtained in the fit 



^^^'l^^'^.t L l^t P^'^'i^'^'^^ °^ '^'^^'^^f''^''^; to the time integrated Dalitz distribution of right-sign 

K~TT^TT^ decays. The relative phase Sj also de- 



eters is hidden in Ai_2 = wi,2 — i^i,2/'^- If CPV 
is neglected the amplitude for tagged decays is 
7V((m^, m?_, t) = (m?_ , , t) . Amplitudes for D 
decays are parametrized in the isobar model as a sum 
of 18 Breit-Wigner resonances and a constant non- 
resonant term. The result of the fit in terms of mixing 
parameters is presented in Fig. [Sj 

Numerically, the fit which allows for the CPV re- 
sults in 

X = (0.80 ± 0.29 ± ^■^^)% 
0.16 

(0.33 ± 0.24 ± ^■!*?)% 



y 
k/p\ 



pends on M^{K~^Tr~), ]VP{K+Tr°), and is determined 
from the fit to the wrong- and right-sign Dalitz distri- 
butions. Parameters x" and y" are, similar as in the 
case of — > K^tt~ decays, a rotated mixing param- 
eters X and y, now by an unknown strong phase shift 
S^^T^o between two points in phase spaces of DCS and 
CF decays to R-n+n^. 

A fit to the time evolution of the wrong-sign Dalitz 
distribution results in 

^7 

x" = (2.61 ± ±0.39)% 



0.67 
55 

y" = (-0.06± • ±0.34)% 
0.64 



(16) 
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2.4. V(3770) ^ D^D^ 

Pairs of neutral D mesons are produced at the 
threshold in a coherent C — ~1 state. With the Cleo-c 
detector one exploits this quantum coherence to deter- 
mine several parameters related to the mixing. The ef- 
fective branching ratios for ^'(3770) D'^D^ /1/2 
are modified with respect to the values measured from 
uncorrelated meson decays. The branching frac- 
tions can be obtained by integration of the time de- 
pendence for V'(3770) ^ D°D° /1/2 [H, 

+ (|a+p - |a_p) cosxAt- 

- 23fi(a;a_)sinhyAt-|-23(a;a_)sinxAt , (17) 

where At is the time interval between two D mesons 
decays, and amplitudes a+ and a_ are in the limit of 
no CPV defined as a+ = Af^Af^ — Af^Af^, a_ = 
Af^Af^ — Af-^Af^. The branching fraction for de- 
cays to the final state composed of a CP-even state 
(S'+) and semileptonic state {e~X), for example, is 
found to be Sr(V'(3770) ^ D"D° S+e~X) w 
Br{D° S+)Br{D° e-X){l - y). In a similar 
manner other double-tagged (both /i and /2 recon- 
structed) branching fractions depend on the parame- 
ters y, Rm — (x'^ + y^)/2 and \JRd cos 5 Single- 
tagged (only a single final state reconstructed) branch- 
ing fraction {Br{D° S+) and Br{D° e- X) in 
the above example) remain unchanged compared to 
the uncorrelated measurements. 

The measurements consist of determination of a set 
of double- and single-tagged branching ratios. Ex- 
amples of reconstructed double-tagged signal yields 
are shown in Fig. [G] Parameters mentioned above 
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Figure 6: Fully reconstructed V(3770) D° D° /1/2 
decays ^15]. The hatched area represents events outside 
the signal region. 

are obtained from a simultaneous fit to the measured 
single- and double-tagged branching ratios (includ- 
ing additional world-average information on individ- 
ual branching fractions). The result is 

y = (-5.2 ±6.0 ±1.7)% 



Rm = (2.0 ± 1.2 ± 1.2) x 10"^ 
^cos5 = 0.089 ±0.036 ±0.009 . 



(18) 



This is the only direct measurement of the strong 
phase difference 6 between DCS and CF decays — ^ 



2.5. Average of results 

Measurements presented in the previous sections as 
well as others constrain the possible values of the 
mixing and CPV parameters in a way specific to the 
decay mode under consideration. A schematic view of 
constraints on x and y posed by some of the measure- 
ments is presented in Fig. [T] Heavy flavour averaging 
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Figure 7: Illustration of constraints in the {x, y) plane im- 
posed by measurements in individual decay modes. Rm 
corresponds to measurements of semileptonic decays. 

group [16| performs a fit to the D° mixing and 
CPV related observables. Correlations among the 
measured quantities are provided by the experiments 
and included in the fit. The average values from the 
fit in which CPV is allowed for are presented in Tab. 

ID 

Table I World average values of mixing and CP violation 
parameters [ig[. 



Parameter 


Value 


X [%] 

y [%] 


0.89 ± If^ 
0.75 ± III 


5n 


21.9 ± II I 


Rd [%] 


0.3348 ± 0.0086 


Ad [%] 


-2.0 ±2.4 




0.87±°:i« 




-9.1 ±«:^ 




no Q 1 25.9 



The no-mixing scenario {x, y) — (0, 0) is excluded 
by the world average of results at the level correspond- 
ing to almost 7 standard deviations. Since positive 
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values of x and y are preferred the almost CP even 
state of neutral D mesons seems to be shorter-lived 
and heavier. 



3. Prospects and summary 

A rich harvest of results in the field of mixing, 
arising from a range of various experiments, only in 
the last two years established the oscillations of this 
neutral mesons with a significance of around seven 
standard deviations. The measured mixing param- 
eters X and y of 0(10"^) are at the upper edge of 
the values that can be accomodated within the SM. 
The results impose severe constraints on the parame- 
ter space of a wide range of New Physics models • 
To illustrate this, Fig. [5] shows the constraints on the 
squark mass - coupling constants plane arising from 
the established value of x in the supersymmetric R- 
parity violating models. 




\u Excluded 



0.02 0.03 0.04 



Figure 8; Constraints in squark mass - R-parity violating 
couplings plane arising from the established value of x. 
The hatched area represents the excluded region from x < 
1.2 %. The lines represent excluded regions for x below 
1.5% - 0.3% 

Several measurements focus on the search for CP 
violation in the D'^ system where a positive signal at 
the current level of sensitivity would represent a clear 
indication of New Physics processes. At the moment 
there is no such hint, and CPV is in several processes 
found to be below the 0.3% level. 

In the near future the major experimental task in 
the field is to measure x with a higher accuracy and to 
further limit the range of the CP violation in various 



decay modes. Recent experimental results might trig- 
ger further theoretical efforts in the SM predictions 
for the mixing parameters, although at the moment 
more precise calculations appear to be difficult. 

A more or less educated guess tells that with a mod- 
est integrated luminosity of 5 ab~^ collected by the fu- 
ture Super-B factory it would be possible to measure x 
and y with a precission of 0.1% - 0.15%. Furthermore, 
the CPV parameters \q/p\ and (j) could be measured 
to around ±0.1, enabling tests of the CP asymmetries 
at around 10"'^ level. 
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